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SUMMARY 


An extended form of the Ackeret iteration method, 
applicable to arbitrary profiles, is utilized to calcu- 
late the compressible flow at high subsonic velocities 
past an elliptic cylinder. The angle of attack to the 
direction of the undisturbed stream is small and the 
circulation is fixed by the Kutta condition at the 
trailing end of the major axis. The expression for the 
lifting force on the elliptic cylinder is derived and 
shows a first-step improvement of the prandtl-Glauert 
rule. It is further shown that the expression for the 
lifting force, although derived specifically for an 
elliptic cylinder, may be extended to arbitrary symmetrical 
prof lies . 


INTRODUCTION 


Two methods of approximation, complementary in 
nature, have been used for the solution of problems of 
flow past bodies in the subsonic-velocity range. They 
are the Janzen-Eayleigh method (reference 1), in which 
the velocity potential or the stream function is developed 
in a power series of the stream Mach number, and the 
Ackeret iteration method (reference 2), in which the 
velocity potential or the stream function is developed 
in a power series of a geometrical parameter character- 
istic of the body. The complementary nature of these 
two methods lies in the fact that the ' Janzen-Rayleigh 
procedure yields accurate results in the case of thick 
bodies, for which the critical stream Mach numbers are 
low; whereas the Ackeret iteration process yields accu- 
rate results in the case of slender shapes, for which 
the critical stream Mach numbers are in the neighborhood 
of unity. 



2 ' 


MCA TN Mo. 1113 


In the Ackeret Iteration process the assumption is 
made that, if X is a parameter that represents the 
departure of the profile shape from a straight-line seg- 
ment at zero angle of attack (that is, X may be the 
thickness, the camber, or the angle of attack), the 
stream function x|/, say, may be developed in a power 
series of X in which the coefficients are functions of 
the plane coordinates x and y and of the stream Mach 
number 1\U . ‘That is, 

-L 7 

^ = -Uy + f-^XjyjMjJx + f 2 (x,y;M 1 JX 2 + ... (1) 

where U is the velocity of the undisturbed flow. By 
inserting this expression into the fundamental nonlinear 
differential equation for \|/ and by equating the coef- 
ficients of the various powers of X to zero, a system 
of linear differential equations for the functions 

f n ix,y;M-^ ; is obtained. The integration of this system 

of differential equations cai be performed for as many 
steps as desired, the first step f^^yjM^j being 

equivalent to the Prandtl-Glauert approximation. There 
is, however, a fundamental difficulty with the form of 
the development, equation (1),. which does not appear in 
the Janzen-Rayleigh method. In the Janzen- Rayleigh 
method the expansion of the stream function in powers of 
the stream Mach number, namely, 

\J/ = \lf Q + ^ 1 M 1 2 + + ... (2) 

can always be obtained, whereas the possibility of the 
expansion, equation (1), cannot be guaranteed beforehand. 
This difficulty in the Ackeret process was avoided in 
references J and i| by choosing as solid boundaries pro- 
files having no stagnation points. In such cases, a 
development of the form of equation (1) is always possible. 
When, however, shapes are chosen that possess stagnation 
points, terms of the type \ n log X ultimately appear on 
the right-hand side of equation (l) and the explicit 
development in powers of "X is strictly no longer pos- 
sible. This difficulty may be avoided by assuming a 
somewhat more general form for the development of 
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the stream function; namely, 


\| / = ~Uy + 'l' 1 (x,y) + ^ 2 ( x >y) + • •• (3) 

where the shape parameter X is contained implicitly in 
the various functions \j/ n . In equation (3) corre- 
sponds to the Prandtl-Glauert approximation and is 

made to satisfy the exact boundary conditions at the 
solid and at infinity. For the purpose of defining and 

controlling the iteration procedure, the function \|/ , „ 

° n+1 

is regarded as small compared with the preceding func- 
tion and the derivatives have a similar relationship. 

It can be stated now that the aforementioned difficulty is 
more apparent than real, for in carrying out the iterative 
steps according to equation (1) X behaves like a dummy 
parameter, which serves only to regulate the iteration 
process in exactly the same manner as equation (3), 

For slender bodies the first few steps of the Ackeret 
iteration method may be expected to yield a good result 
with the exception of a small region in the neighborhood 
of a stagnation point. Hantzsche, in reference 5 > treated 
the case of an elliptic cylinder at zero incidence in a 
uniform stream according to Ackeret* s process as repre- 
sented by equation (1), with the exact boundary condi- 
tions at the solid and at infinity being satisfied. 

These calculations were carried through the X3 

-terms 

and included a term X^+ log X. A comparison, made in 
reference 5 > of this result with that obtained by the 
present author according to the Janzen-Rayleigh method 
showed agreement in the terms common to the two develop- 
ments. A similar comparison made for the "bump" in 
reference 5 also showed complete agreement in the terms 
common to the two methods. These comparisons illustrate 
the fact that solutions obtained by the Janzen-Rayleigh 
and the Ackeret methods are simply different representa- 
tions of a single unique solution. In particular, it is 
evident that the Janzen-Rayleigh method, which is valid 
at stagnation points, includes parts of the Ackeret 
development and, conversely, the Ackeret development 
includes parts of the Janzen-Rayleigh result. Although 
the first f ew terms of equation (3) do not yield very 
accurate results at stagnation points, these terms 
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nevertheless represent correctly, to some extent, the 
effect of compressibility at such points. The accuracy 
of the calculations obviously depends on the number of 
terms \[ f n derived, each new term reducing the extent of 

the region of inaccuracy in the neighborhood of a stag- 
nation point. 

The question of the convergence of the sequence of 
functions, equation ( 3 ), is a. difficult one and should 
be thoroughly investigated. Schmieden and Kawalki 
(reference 6; state that both the Janzen-Rayleigh and 
the Ackeret developments diverge when the local velocity 
of sound is first exceeded in the region of flow, that 
is, at the critical value of the stream Mach number. 

Ttiis statement is contradicted, however, by results 
obtained by means of hodograph or velocity variables. 

In general, as shown in reference 7> the limit of steady 
potential flow of a compressible fluid is determined by 
the vanishing of the Jacobian of the transformation from, 
physical x,y- variables to hodograph 9 ,q- variables . The 
vanishing of this functional determinant does not neces- 
sarily occur when the local Mach number first reaches 
unity and, consequently, there exist continuous solutions 
of the general differential equation governing the flow 
of a compressible fluid for which a part of the region 
of flow is supersonic. It is reasonable to conclude, 
therefore, that the series solutions given by equa- 
tions (2) and ( 3 ) diverge at the value of the stream Mach 
number that marks the limit of potential flow rather 
than at the value for which the local velocity of sound 
is first exceeded in the flow. 

The Ackeret process in the form of equation (1), but 
for the velocity potential, was used to calculate the flow 
past a bump and past a circular arc profile in refer- 
ences 3 f *nd 4, respectively. This calculation was 
possible because the profiles considered did not possess 
stagnation points. As a consequence, the problem* could 
be treated completely in the physical .plane and, more- 
over, the boundary conditions on the velocity potential 
were tractable. In general, however, ’when shapes with 
stagnation points are treated, the Ackeret process in the 
form of equation ( 3 ) instead of in the form of equation (1 
should be utilized. In such cases the problem is most 
conveniently treated in a new plane related to the plane 
of flow by an affine transformation. In this plane, 
however, the boundary conditions for the velocity poten- 
tial become very complicated. On the other hand, the 
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boundary condition for the stream function, \|/ = 0 at 
the solid, is invariant for affine transformations, and 
it is therefore suggested that the stream function be 
utilised. The transference of the problem to an affinely 
connected plane unfortunately introduces a distortion in 
the actual profile that varies with the stream Mach num- 
ber. In the case of an elliptic profile, however, this 
distortion does not matter, for the affine distortion of 
an ellipse leads to another ellipse. This property of 
the ellipse makes it a preferred profile for the Ackeret 
iteration process in the form of equation (5) and is also 
the reason it is chosen as the solid boundary for the 
problem treated in the present paper. Specifically, the 
problem treated herein is the determination of the flow 
of a compressible fluid past an elliptic cylinder set at 
a small angle of incidence in a uniform stream, with 
circulation determined by the Kutta condition at the 
trailing end of the major axis. The main purpose of this 
calculation is to obtain some Information with regard to 
the effect of compressibility at high subsonic stream 
Mach numbers on the lifting force acting on an elliptic 
cylinder. A calculation is now in progress in the 
Physical Research Division of the Langley Memorial Aero- 
nautical Laboratory to find the effect of compressibility 
at high subsonic stream Mach numbers on the moment and 
on the location of the center of pressure of an elliptic 
cylinder. 


CALCULATION OP THE FIRST AND SECOND APPROXIMATIONS 


The equation of continuity can be written as 



and the condition for irrotational flow as 


dv bu _ 
dX ~ dY 


( 5 ) 
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where 

X, Y rectangular Cartesian coordinates in physical flow 
plane 

u, v components of velocity along X-axis and Y-axis, 
respectively 

p variable density of fluid 

constant density of undisturbed fluid at infinity 
Equation (5) defines a velocity potential 0, where 


/ 






(6) 


and equation (Iq) defines a mass flow or stream func- 
tion if , where 


— — Pq d\j/ 

" 5x “ p 5 y 

> (7) 

_ 60 _ Pq 6\1; 

" bY ~ “ p 6X 


From equations ()q) and (5), with the use of equations (6) 
and (7), respectively, the following alternate forms of 
the basic differential equation are obtained: 


q_/£_dA + W_e_6A = 0 

oX \ 9 1 6x ) 8y \Pq 6y ) 

and 


6 /Pq 6\j/ \ + 6 / P 1 d\jj \ _ q 
6x \ P 6xy 6 y \ P 6y/ 


( 9 ) 
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For reasons stated in the "Introduction," equation (9) 
for the stream function will be treated in the present 
paper. The difficulty of obtaining a solution lies in 
the fact that the density of the fluid p is related 
to u and v because of the dependence of p on the 
local pressure. If it is postulated that the fluid is 
nonviscous and if the fluid is assumed to flow from 
infinity with a constant velocity U, the total energy 
of the fluid will then have the same value at every point 
in the region of flow. If p and q denote the pres- 
sure and the velocity of the fluid, respectively, 
Bernoulli’s equation becomes 

SE + i(q 2 - U 2 ) = 0 (10) 



where the lower limit of integration refers to the 
starting conditions at infinity. Moreover, p and p 
are connected by some known adiabatic equation of state 
such as 




Constant 


(ID 


where, in the case of a perfect gas, y denotes the 
ratio of specific heats at constant pressure and at 
constant volume. 

By means of the definition of the velocity of sound. 



( 12 ) 


the following relations are obtained from equation (10): 
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P 




(15) 


where M, = — is the Mach number of the undisturbed 
1 C 1 

^ p p 2 

stream at infinity and q- = u + v . 

By means of equations (7 ) > equation ( lit) can t> e 
rewritten in the form 







(16) 


where \|/ x and \|/y denote bty/bX 
tiveiy. If, for the moment, Pq/P 


and 

bty/bY, 

respec- 

and. 

<; x 2 + •• 

it. i 

T7 2 



are considered to be dependent and independent variables, 
respectively, a Maclaurin expansion in the neighborhood 


of the undisturbed stream, where 

Pq 

— = 1, yields 
P 


+ MV? 


X- -1 



0 and 
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of this development has been chosen to be consistent 
with the Ackeret iteration process, which is essentially 
an iteration around an undisturbed stream of zero inci- 
dence. Corresponding to equation (3), 

\J/ = -UY + \|/^ + ^2 + • • • (18) 

where ^ n+ ]_ is regarded as small compared with and 

the derivatives have a similar relationship. When this 
expression for ^ is substituted into equation (17) and 

2 2 

■'I'p-v- is noted to be of the same order as or , 

^ IX 1Y 

then inclusive of the second power in the derivatives, 

T ■ 1 - - OV • (* 2 ' *)(¥ • 

2 

+ |([i 2 - l) 2 |"(Y + k) + (Y + 1)(m - 2 ” l) + • • • (19) 



When this expression for P^/P - s substituted into the 

basic differential equation (9) and terms of equal order 
of magnitude in the derivatives of are equated, the 

following differential equations for and i| are 

obtained: 
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These differential equations are most easily solved by 
first applying the affine transformation 


A 


x = X 


y - i-Y = Jl - M-, 2 Y 

\X V X 


> 


Equation (20) then becomes a Laplace equation 


^lxx + ^lyy ® 


and equation (21) becomes a Poisson equation 

r 


if 


2 xx 


^2yy = (p 2 - i) \ 


~ ^ v u** +2 y7%r 


(f - 0 

1 

pu !yy 


J 


( 22 ) 


(25) 


(2^) 


For purposes of calculation a new stream function ^ 
may be conveniently introduced, where 


^ = pi#* 

From equations (22) it can be seen that 

■ *X = *x = 


(25) 


tyy = i \J/ = UVj/ 


n y 


and 
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so that the undisturbed flows, at zero angles of inci- 
dence, in the physical and affinely distorted planes are 

. With the intro- 

1 

duction of the stream function equations ( 23 ) 

and ( ) become 


identical; that is 


(h), = O 


+ \J/* - 0 

ixx y lyy 


( 26 ) 


and 


r * 2 xx + ** 2 yy = { *\y**lxx + ^'"lx^'lxy 


5 + 


( Y + l)(p 2 - l) 


^ iy^ lyy 


1 


or, with the use of equation ( 26 ), 


'I'W'i'hyy 




2 + ( Y + 1 ) (p 2 - l) 


^ ly^ lxx 


+ ixy 


(27) 


Equations (26) and ( 27 ) show that, in order to cal- 
culate the various approximations, the incompressible 
flow past the distorted profile in the xy-plane must be 
known. Thus, if a profile is given in the physical 
XY~plane, it will be necessary to find the conformal 
transformation to a circle of the distorted profile in 
the xy-plane for each value of the stream Mach number, 
since the affine distortion of a profile depends on the 
stream Mach number, as shown by equation (22). In 
general, then, the problem to be solved is the flow past 
an arbitrary profile in the affinely distorted xy-plane.. 
This procedure in general involves the laborious calcu- 
lations of the coefficients of conformal transformations 
to a circle for a number of values of the stream Mach 
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number in such a way that the distorted profiles in the 
xy-plane correspond to the given profile in the physical 
XY-plane. In the case of an ellipse, however, the dis- 
torted profile is again an ellipse, and it is therefore 
a simple matter to transfer the results obtained in the 
xy-plane to the XY-plane of the original elliptic 
profile. In the present paper the elliptic profile is 
so oriented that its major axis lies along the X-axis. 
The relation between the profiles in the two planes is 
then given by 


a = a' 


b = lb' 

H- 

? .2 , - l-u.2 

^ — b» 


R = 


ri + LzJV 

2g 


> ( 28 ) 


tan a = 1 tan a> (or, to the first order, a = - a' ) 


where 


r = 1 (a + b) 


semimajor axis of ellipse in affinely distorted 
plane 


semiminor axis of ellipse in affinely distorted 
olane 


c semifocal distance of ellipse in affinely distorted 
plane 

a angle of incidence of ellipse in affinely distorted 
plane 

and the prime indicates corresponding values in the 
physical plane. 
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The solution of equation (26) is in general the 
imaginary part of an analytic function w-^(z) where 

z = x + iy. (The asterisk has been dropped. ) It is 
easy then to verify that 

^x = ^y( w lz " ? lz) 

*iy = a( w iz + v7 i?) 

V lxy “ 2 "C W lzz + w izr) 

^lxx ~ ~^lyy ~ 2 i( w lzz “ w lzz) 


where a bar indicates conjugate-complex quantities. 
Therefore , 


^l/lxx = T - ? -0 


,) = (« 


w lzz w lz + W lzz w lz) = ( W T- + w t-H*P‘ w V 


lz v 'lz 


a 'lx'li X y =I,P ‘(’ 


*12Z V '1Z -'''lZZ W 


? lz) = - 1 ( w lz - w lz) R - P - w 


Equation ( 27 ) can then be written as 


.2 

o ^2 


K W lz 2 ) + + ^( w lz*lz) (29) 

dzdz ° L z ' z 

where a = (y + 1 ) - l) . 


Consider now a circle of radius R, v/ith its 
center at the origin of the z» -plane, into which the 
distorted profile in the z-plane is mapped by means of 
a conformal transformation. Any point on the circle 
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can be expressed in the form 

z» = x» + iy ' = R( cos £ - i sin £) = Re“ i ^’ 

so that the point describes the circle in the clockwise 
sense and leaves the region outside the circle on the 
left. If now £ = £ + it), the transformation 


z' = Re** 5 -^ (50) 


yields a circle of radius R when r) = 0 • and the 
infinite region of the z' -plane when t) - + 00 . Equa- 
tion (30) can be looked upon as the transformation from 

Cartesian coordinates . x'.,y' to polar coordinates Re 1 !, -£. 
The conformal transformation of an arbitrary profile in 
the z-plane into a circle of radius R in the z» -plane, 
with center at the origin, can therefore be written 

2 = f^e 1 ^ (31) 


where ^ = 0 yields the parametric equations of the 
profile in the z-plane. 

When t and I are introduced as independent 
variables, equation (29) becomes 



(32) 


The general solution of equation (32) or of equation (29) 
is obtained directly and is 


*2 = -j^ 2 -i) I - p - 


2<Jz w lz 2 t (a+]q)w 1 w lz + F(£) 


(33) 


where P(t) is an arbitrary function to be determined 
according to the boundary conditions at the solid surface 
and at infinity. 
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Consider now an ellipse in the affine z-plane with 
semimajor axis a and semiminor axis b. The undis- 
turbed stream at infinity makes a small angle a with 
the negative direction of the x-axis, and the circulation 
is determined according to the Kutta condition that the 
dovmstream end of the major axis is a stagnation point. 

In accordance with the Ackeret process, 

^ = -y + ^ + \j/ 2 + ... (34) 


so that the stream function is expanded around an undis- 
turbed stream at zero incidence. Since the angle of 
attack has been assumed small, powers of a higher than 
the first are neglected. The boundary conditions are 
then, at the surface of the ellipse, rj = 0, 


*1 = y I 

->2 = 0 J 

and, at infinity, r) = 03 > 

^1 _ 

bx 


and 



bx 


0 


> 


H 2 


(35a) 


(35b) 


Now, the transformation 
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maps the region external to the ellipse in the z-plane 

into the region external to a circle of radius R = i.(a+b) 

with center at the origin of the z» -plane. The complex 
potential for the flow considered is given by 

w = «/zte ia + 

\ z 


- - 


2iR sin 


-t z t 
a log _ 


(37) 


When the variable £ is introduced by means of equa- 
tion (JO), equation (56) becomes 


z = c cos (£ + iX) 
and equation (57) becomes 


(58a) 


w = -2R cos (£ - a) - 2 r£ sin a 


(58b) 


where a = c cosh X and b = c slnh X. This expression 
for w includes the uniform undisturbed stream -z which 
must be extracted in order to obtain Thus, for a 

small angle of attack a. 


W]_ = c cos (£ + iX) - 2R cos £ - 2R (sin £ + £)a (59) 

and 

\ = I.P.w-l 

= -c sinh (r) + X) sin £+2R sinh sin £ 

r- 2R (sinh r) cos £ + r])a 

It can be easily verified that this expression for 

satisfies the boundary conditions stated in equations (55) 
and also that the downstream end of the major axis is a 
stagnation point. 

Similarly, the most direct way to determine ^ is 
to consider it to be the imaginary part of a nonanalytic 
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function W£ of £ and £ . 
( 38 ), and ( 39 ), 



- 2— cos £ - 2— a (sin £ 
c c 


Thus, from equations 



(35), 


(I - iX) 


In order to satisfy the boundary condition ^2 = 0 
the surface, rj = 0, it is a simple matter to supply 
the functions of £ needed to make the coefficients 
of w and w lz 2 vanish for r) = 0. For example, 

cos (£ - iX) = cos (£ - iX) for n = Thus, 



+ ( o + l+JcTcos ( £ - iX) - 2— cos £ - 2— a (sin £ + £) 
I c c 


-cos ( £ - iX) + 22. cos £ + 22a (sin £+ £) 


w lz > (1+0) 


In order to satisfy the boundary condition at infinity 
and the condition that the downstream end of the major 
axis be a stagnation point, the procedure, according to 
equations (7), (22), and (25 ), is as follows: 


Pl U 




21 d£ 6 ^2 _ d£ 6 
dz 6£ dz 6 £ 




( 1 * 1 ) 


From equations (39) to (1+1) it then follows that at 
infinity 
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JL_ 

P x u 






+ « rb ia 


Hence, in order to satisfy the boundary conditions at 

^ o „ 6\(/o 

infinity, u = —• = 0 and v 0 = £. = 0, a term 

i r ? \ 7 6x . 

“I7 b v “ -v( c + s ^ n ^ imaginary part of which 


vanishes for r) = 0, must be added to the expression 
for w 2 . This addition to w 2 introduces a velocity at 
the downstream end of the major axis, given by 




£=tt, T)=0 


-Ki 


„ 2 - 



Again, in order to render (£ = tt, t} = 0) a stagnation 
point, a term -jjb(ji^ - l)(P + l|)a£ satisfying the 

boundary conditions at the solid and at infinity must be 
added to the expression for w 2 , Finally then 


( l - iX ) 

(I - iX ) «2~ cos I 

- 2 ^. 0 . (sin £ + £ ) - cos (£ - iX ) +2 — COS £ 
c o 

+ 2^-a (sin s + n| w i z + 2b(a + (sin £ + £ ja j (Ij2) 



The complete expression for \|/, obtained from equa- 
tions (3l^), (39 ), and (k2), is then given by 
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2R cos £ + 2Ra ( sin £ + £ ) 


+ |(p, 2 - l)^°c jcos ( £ - iX) - cos ( £ - i\)J v/ lz 2 
+ (a + i|)c£cos (£ - iX) - 2~ cos £ ~2~a (sin £ + £ ) 


- cos (£ - iX) +2^- cos £ + 2-a ( sin £ + £ )J 


w 


lz 


+ 2b (a + 1|) a (sin £ + £ ) 


where w-^ is obtained from equation (39) and 


( 1 * 5 ) 


w lz = 


dw x d£ 

d £ dz 


CALCULATION OP TIIE VELOCITY COMPONENTS u AND v 


The components u and v of the velocity of the 
compressible fluid past the actual profile in the physical 
XY-plane can be put into the following convenient form by 
means of equations (7)> (22), and ( 25 ): 




dx 


= 21 ai = at i_ ( w . w> Ui4) 

dz d£ dz bi 


where w is the expression on the right-hand side of 
equation (Ijy ) with I.P. omitted. It follows that 
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where equations (28) provide the correspondence be tween_ the distorted 
ellipse in the z -plane and the actual ellipse in the Z -plane . 
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Equation (45) will no w be utilized to calculate the 
lift on an elliptic cylinder in compressible flow. For 
this calculation a control contour, which eventually is 
taken as infinitely large, will be applied in the usual 
manner. The elimination of the variable £ in equa- 
tion (45 ) & nd in the equations needed for the calculation 
of the lift is in general impossible and in the present 
example undesirable. Since the regions at infinity corre- 
spond in the z- and z> -planes, it is convenient to 
choose a large circle in the z* -plane as the control 
contour and to effect all the calculations with £, and r) 
as independent variables. The advantage gained by this 
procedure is a great reduction in labor, in that on a 
circle. r) is constant and hence only functions of the 

( 

single variable z* e J appear. The first step is 
to obtain developments for pu/p-^U and pv/p-^U in the 

neighborhood of infinity and then to form the combina- 
tion -2— (u - iv). This calculation, according to equa- 

p-jU 

tion (45), yields the following result! 


(u - iv )=-(!+ i pa ) + 


pi D 


xa 

a *1 


+ p(t- 

2 -l)(a + l±) b 

4 V 

y ' a+ b 


L. 




i~_A_ 

j|_a + b 


+ K ^) 2 ^- 1 ) (0 + w 


(1 - p)z t 2 + (1 + p) ^ 


+ ia 


- xa 


cl 


a + b 
b 

a + b 


- (1 + p) 


(>' - ■) 

hi ; 


2 


(1 - p)z t 1 - - (1 + p) 


= -2J 
1 


p/zt 2 + + (1 - p)z A 

Z, V 


'MJ 


• • • 


(46) 
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where z» = e -1 ^. Since and 

it follows from equations (1?) an d (i+6) that 



CALCULATION OF THE LIFT 

In a compressible flow as in an incompressible flow, 
the lift is given by (reference 8) 

l c = Pl nr 0 (V3) 

where T is the circulation round the profile and 
c 

where, by definition, 

r =n(u dX + V dY) = R.p. Iff) (u - iv)dZ (4-9) 
c J o 

Nov/, from equation (22), 

X = x 

y = . wr 
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Hence 


Z = L±Ji z + L-. Fa 
2 2 

Since z = c cos (£ + iM> 

dZ = sin (£ + iX)d£ - -^-^ c sin (£ - i\')dl 


For the evaluation of the line integral for T , the 

C/ 

control contour is a large circle in the s' -plane. 
Therefore, n is a constant and d£ = d £ = i~; hence. 


dZ = - 


ce 


T]+\ 


(1 + 1 - 0 / 


*»e 2r ' +2X 


- z t 


+ (1 - [i)( 


z t 


Z 1 




dz 1 

z' 


(50) 


— j_£ 

where z» = e . The desired expression for u - iv, 
obtained from equations (J.|.6) and ( iq.7 ) , is 


— (u - iv) = -1 *- i(xa+ — 
U e r) 


1 + p(|j 2 - l) (a + i|.) — 
Ip ' a + 


(1 - h) z ’ 




+ • • . 


(5D 


Equations (I 4 . 9 ) to ( 51 ), with only terms containing the 1 
factor dz*/z» contributing to the line integral for T c , 
then yield the following result: 
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2l+ 


r = [(.irRU [i a 


1 + K“ 2 - O* 0 + 


(52) 


If a, b, R, and a are replaced by a», b», R», and a» 
according to the correspondence equations ( 28 ), then for the 
actual elliptic cylinder in the physical flow plane 2f, 


r 


c 


lpnR T Ua» ix + 2irUa»b* 





Since the circulation in the incompressible case, 
M, =0 or p = 1, is 


r i - IlttR * Ua* 


the ratio F c /r^ is given by 


i 2R* 


(1 - p) + r(| m, 2 - l)(o‘+ k) 


With R» = i(a» + b») and a = ( Y + l)(p 2 - l), 


^c = ^c = . 

■r T*1 ^ 


t « 


r. 


i + t« 


m-(m- “ 1 ) + ^(y + 1 ) (m 2 - l) 


(55) 


where t’ is the thickness coefficient b’/a* of the 
elliptic cylinder in the physical flow plane. This equa- 
tion represents a first-step improvement of the Prandtl- 
Glauert approximation and reduces to that result when 
t»-->0. 

Although equation (53) has been derived specifically 
for an elliptic cylinder, it will be shown that the result 
can be extended to a slender arbitrary symmetrical profile. 
Hantzsche and Wendt (reference 9) derived a similar 
relation for the case of a symmetrical Joukowski profile 
with a sharp trailing edge. The result of Hantzsche and 
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Wendt may be written as 



y. + 


1 + € 


[i([i - 1) + i( Y + l)(p 2 



(5k) 


Note that the function of jjl contained between the 
brackets is the same in equations (53) and (5M* This 
coincidence suggests a correspondence between the 

factors — — and — € . 

1 + t» 1 + e 

A correspondence is obtained in the following 
manner: It is well-known that by means of the mapping 

function 

c 2 

Z = Z» + — 

li/zt 


the circle of radius c^/k > with its center at the 
origin of the Z’ -plane, is mapped into the line segment 
extending from s = -c to s = c, and the circle of 

radius —(1 + €'), with its center at Z* = e— , is 
2 2 

mapped into a symmetrical Joukowski profile with sharp 

trailing edge in the Z* -plane. Now 

( 55 ) 

1 + t* 2R 


or, with R = 


¥ 1 + c) ’ 


1 + € - 


t* 


1 + e 


2(1 + €) 

= _J_ . I ( '_£_V 

1 + € 2 \1 + v 


1 + t* 
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4 - i € 

Thus, to the first power in € , — - = — and 

1 + t» 1 + c 

the correspondence between equations (53) and (54) is 
established. 

In the case of an arbitrary symmetrical profile, 
the Theodorsen method (see reference 10) is particularly 
well suited to obtain an expression corresponding to 
t * ' 

. An essential feature of the potential theory 

1 + t* ■ 

of arbitrary wing profiles developed by Theodorsen is a 
rapidly convergent procedure for obtaining the conformal 
transf ormation of the profile to a circle, also the 
radius R of the circle. The coefficient of the 

— — term of this conformal transformation, denoted by 

Z* 

c 2 /k, and the radius R of the conformal circle define 
an ellipse 

Z = c cos (£ + i^< ) 

with 


Then, from equation (55 )> 


1 1 

1 + t» 




and, therefore, for an arbitrary symmetrical profile, 
the formula that corresponds to equation (53) may be 
written 






(56) 


Table I shows values of the ratio L C/ /Lj_, calcu- 
lated by means of equation (53 )> f or> various values of 
the thickness coefficient t* and the stream Mach 
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number (with y = 1*4 for air) and figure 1 shows 

the corresponding graphs with M-j_ as abscissa and 
L^/L i as ordinate. 

Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., May l6, I 9 J +6 
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